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ABSTRACT 

A follow-up survey using the Submillimetre High-Angular Resolution Camera 
(SHARC-II) at 350 /im has been carried out to map the regions around several 850 /im- 
selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). 
These observations probe the infrared luminosities and hence star-formation rates 
in the largest existing, most robust sample of submillimetre galaxies (SMGs). We 
measure 350 /im flux densities for 24 850 fim sources, seven of which are detected at 
> 2.5cr within a lOarcsec search radius of the 850 /im positions. When results from 
the literature are included the total number of 350 /mi flux density constraints of 
SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far- 
infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich 
(Mdust — 9 x 1O 8 M0), luminous (£fir — 2 x 10 12 L q ) star-forming galaxies with 
intrinsic dust temperatures of ~ 35 K and star-formation rates of ~ 4OOM yr~ 1 . We 
have measured the temperature distribution of SMGs and find that the underlying 
distribution is slightly broader than implied by the error bars, and that most SMGs 
are at 28 K with a few hotter. We also place new constraints on the 350 /im source 
counts, iV 35 o(> 25mJy) - 200-500 deg~ 2 . 

Key words: submillimetre - surveys - cosmology: observations - galaxies: high- 
redshift - galaxies: formation - galaxies: starburst 
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1 INTRODUCTION 



The SCUBA HAlf Degree Extragalact ic Survey (SHADES; 
iMortier et"al]l2005l : [Coppin et al.ll2006h mapped ~ 0.25 deg 2 
of sky with an RMS of 2 mjy at 850 fim with t he Submillime- 
tre C ommon-User Bolometer Array (SCUBA; iHolland et al.l 
1 19991 ). The area was split approximately evenly between 
the Lockman Hole (LH) and the Subaru- XMM Deep Field 
(SXDF). Using uniform selection criteria, the survey un- 
covered 120 SMGs with a me dian deboosted flux density 
of ~5mJy (|Coppin et al.ll2006l ). The SHADES programme 
was designed to study the nature and evolution of high 
star-formation rate (SFR) submillimetre galaxies (SMGs) 
via a systematic study of a well-characterised and statisti- 
cally meaningful sample. The programme includes an effort 
to identify members of the source list at other wavelengths 
using deep follow-up data from the radio Ijlvison et al.ll2007l ) 
to X-ray, in order to characterise the SHADES population 
and to probe the variation in the star-formation and cluster- 
ing with redshift. The relatively precise positions available 
from the radio data greatly aid in identifying secure counter- 
parts at other wavelengths, which can then be used to pro- 
vide spectroscopic or photometric redshifts (|Aretxaga et al.l 
2007) and to categorise the sources. 

Even combined with knowledge of source redshift, 
SCUBA 850 /im and millimetre wavelength fluxes do not 
constrain the total dust mass of a galaxy because there is 
an ambiguity between column density and source temper- 
ature. The submillimetre (submm) spectral energy distri- 
bution (SED) of a luminous dusty galaxy arises from the 
re-emission at far-infrared (FIR) wavelengths of absorbed 
optical/ UV radiation from regions of intense star-formation 
(see e.g. ISanders fc Mirabell Il996l and references therein). 
Typically t he dust temperat ure is within a factor of 2 of 
T d ~ 40 K i|Blain et alj l2002h . so the restframe SED peaks 
in the range 60-120 and the SED is almost a simple 
power law at the SCUBA wavelengths and longer. 

At a redshift o f (z) ~2 or 3, typical of SMGs 
jChapman et al.l 120051 ). the peak of the SED is shifted to 
be near 350 /im, and the Su bmillimetre High A ngular Res- 
olution Camera (SHARC-II; iDowell et aLll2003l ) at the Cal- 
tech Submillimetre Observatory (CSO) is very well situated 
to provide the photometry needed to constrain the temper- 
atures, and therefore the luminosities and masses, of the 
SHADES sources. 

We have therefore mapped a subset of the SHADES 
catalogue with SHARC-II. In this paper we constrain the 
FIR SEDs of our sample by fitting modified blackbody 
curves to the FIR photometry of 31 SHADES SMGs includ- 
ing SHARC-II data at 350 ^m. Sections [2] and [3] describe 
the observations and data reduction. Section [4] presents the 
350 ;Um flux densities of the SHADES galaxies and the FIR 
SEDs. Section [5] provides a discussion of the results. Con- 
clusions and future prospects are discussed in Section [6] 
We ad opt the Cosmologica l parameters from the WMAP 
fits in ISpergel et ail (|2003l ): fl A = 0.73, On = 0.27, and 
H = 71kms" 1 Mpc" 1 . 



2 SURVEY DESIGN AND OBSERVATIONS 

2.1 Instrument Description 

SHARC-II is a background-limited 350 and 450 /im 
common-user continuum camera with a 3 arcmin x 
1.5arc min field of view (F OV) at the 10.4-metre CSO in 
Hawaii (|Dowell et al.|[2003r ) . The dish has a very low surface 
error (10.4 /im RMS at 350 um) due to the active Dish Sur- 
face Optimization System (|Leonsl 12005) which corrects the 
primary for surface imperfections and gravitational defor- 
mations as a function of elevation angle during observations 
to improve the telescope efficiency and pointing. The result 
is that the CSO is probably the best telescope in the world at 
shorter submm wavelengths. The resulting beam size (with 
good focus and pointing) is 9 arcsec FWHM at 350 /im. 

2.2 Sample Selection 

120 s ubmm sources have been identified in the SHADES 
fields (|Coppin et al.ll2006l ). Obtaining useful photometry for 
the full set would require several hundred hours of excel- 
lent weather, so we have chosen to observe a carefully cho- 
sen subset of the SHADES catalogue. Many sources lie close 
enough together that multiple targets can be selected within 
a single SHARC-II FOV. We chose eight of our 14 fields 
to contain a large fraction of the SHADES close angular 
pairs. This is mostly just for efficiency, but there is also 
a small chance that the observations would help measure 
clustering in the SHADES catalogue. We chose a deliber- 
ate mixture o f sources with com pact, extended, or no radio 
counterparts (jlvison et al.ll2007l ): 11 sources with one robust 
compact radio ID; three sources with two robust compact 
radio counterparts; three sources with one extended-looking 
radio counterpart; and seven sources with no reliable radio 
ID. Choosing targets with this mixture of radio character- 
istics could help refine the SHADES redshift distribution 
and test if there is a sub-population w ith a high-redshift tail 
fe.g. lDunlop]|200ll . llvison et al.ll2006l) . 

Two of our fields wer e also o bserved by iKovacs et al.l 

|2006t ) or lLaurent et all (|2006l ) (sources LOCK850.3, 
LOCK850.1, and LOCK850.41), which allows for cross- 
calibration and checks of systematics in observation and 
analysis. One field (LOCK44/45) covered two 850 fim source 
candidates from a preliminary SHADES catalogue which 
were later eliminated from the official SHADES cata- 
logue since they we re deemed likely to be spurious sources 
l|Coppin et alj 12006). This field is not used in the ensuing 
analysis. 

2.3 Observations 

Ten SHADES Lockman Hole (LH) fields were mapped over 
the course of four superb stable weather nights in March 
2005 and February 2006 (0.035 < r 22 5 GHz < 0.06). Four addi- 
tional fields were observed in more marginal observing condi- 
tions. For each field data were collected in 10 minute scans 
using a non-connecting Lissajous scan pattern with small 
amplitudes of typically 20 arcsec in both altitude and az- 
imuth and a 'period' of 15-20 seconds. Integration times, 
typical weather conditions and resulting map depths of each 
SHARC-II field are given in Table □ 
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Table 1. Summary of the SHARC-II observations and map properties. RA and Dec. (J2000) coordinates are given for the field pointing 
centre, with each field containing one or more SHADES sources. The date of observations, weather conditions (mean atmospheric opacity 
at 350 /tin), and the 'effective' integration times (see Appendix A; note that this is not the same as the actual time spent on the sky) 
are listed. Since a large fraction of each map is comprised of noisy edges, we quote a median value of a 'trimmed' noise map containing 
only pixels with more than 60 per cent of the maximum depth coverage achieved by the central pixel in the map. Total map areas are 
also given. 



Field 


SHADES 


RA 


Dec 






UT Date 


Mean 


Eft. Int. Time 


Median RMS 


Map Area 


Name 


Sources 


(J2000) 


(J2000) 




(yy-mm-dd) 


r 350 pin 


(min.) 


(mjy) 


(arcmin 2 ) 


LOCK21/28 


21, 28 


10 h 52 to 


57?3 


57° 


30' 


54/.' 


o 


05-02-28/-03-08.09 


1.30 


8.5 


13 


3.1 


LOCK26/32 


26 


10 h 52 m 


39? 7 


57° 


23' 


10': 


'4 


05-02-28 


0.99 


9.3 


9 


3.1 


LOCK4/69 


4 


10 h 52 m 


05? 5 


57° 


27' 


06': 


'4 


05-02-28 


0.88 


9.9 


10 


3.1 


LOCK3/47 


3, 47 


10 h 52 m 


37? 1 


57° 


24' 


57': 


'5 


05-02-28 


0.88 


1.1 


20 


3.0 


LOCK33/42 


33, 77 


10 h 51 m 


56? 5 


57° 


23' 


05': 


'5 


05-03-01 


0.90 


8.9 


9 


3.1 


LOCK10/48/64 


10, 48, 64 


10 h 52 m 


52 s 


57° 


32' 


51': 


'4 


05-03-01 


0.89 


9.8 


13 


3.7 


LOCK1/41/63 


1, 6, 41, 63 


10 h 51 m 


57 s 9 


57° 


24' 


60': 


'6 


05-03-01 


0.80 


13.7 


23 


3.9 


LOCK16/50 


16 


10 h 51 m 


49? 3 


57° 


26' 


44': 


'5 


05-03-08 


1.33 


7.5 


17 


3.8 


LOCK22/25 


22 


10 h 51 m 


35? 2 


57° 


33' 


26'.' 


'9 


05-03-08,09 


1.52 


6.0 


17 


3.5 


LOCK44/45 


none 


10 h 51 m 


56? 8 


57° 


28' 


37': 


'9 


05-03-09 


1.33 


4.7 


22 


3.8 


LOCK15 


15 


10 h 53 m 


19? 1 


57° 


21' 


io': 


'5 


06-02-24 


1.52 


1.1 


55 


2.8 


SXDF1/11 


1, 11 


02 h 17 m 


27?9 


-04° 


59' 


38': 





04-08-30, 31/-09-01 


1.83 


6.4 


14 


3.8 


SXDF3/8 


3, 8 


02 h 17 m 


43? 2 


-04° 


56' 


12': 


'5 


04-08-30, 31/-09-01 


1.81 


5.2 


15 


4.3 


SXDF17/26 


17, 119 


02 h 17 m 


55? 8 


-04° 


52' 


50': 


'4 


05-03-01,03,05,07,08,09 


1.39 


6.8 


17 


3.1 



By design, our strategy was to reach a map RMS of 
3<t ~ 30mJy at 350 /im, sufficient to detect or place useful 
limits on an S S 50 ~ 8 mjy SMG with T d ~ 40 K at z < 3. 
Good weather is scarce, and it is important to be as ef- 
ficient as possible by only integrating down to the planned 
noise level. We tracked the effective integrating time on each 
source by compensating each 10 minute file for zenith an- 
gle and atmospheric opacity at 350 /im, as inferred from the 
7350 and 7225 GHz dippers. See Appendix lAl for details. 
Typically, 3-4 hours of observations were required for a map 
3 arcmin x 1.5 arcmin in size. 

Pointing, focus checks, and calibration were performed 
hourly on standard sources in close proximity to the science 
targets. The same scan pattern was used as for the science 
targets, but with integration times of only 120-160 seconds 
(as typical flux densities are > 2 Jy). 

The typical pointing RMS of 2-3 arcsec at the CSO is 
a non-negligible fraction of the 9 arcsec beamsize, and much 
care is taken to minimise pointing errors. Observations of 
point-like galaxies, quasars, protostellar sources, HII regions, 
and evolved stars are used by the CSO staff for constructing 
a pointing model. The model predictions for the calibrators 
observed before and after science observations are compared 
to the actual pointing measurements; offsets are calculated 
and applied to the model pointing predictions for the sci- 
ence observations during the map coaddition stage of the 
reduction. This procedure yields a pointing accuracy RMS 
of ~ 2 arcsec (D. Dowell, private communication). 

Flux calibration is performed by comparing the known 
and measured flux densities and beamsizes obtained for 
standard calibration sources. For sources in the SXDF, oCeti 
(a Mira variable) and occasionally OH231 (a proto-planetary 
nebula) were used as pointing and calibration sources. For 
sources in the LH, CIT6 (an evolved star) was used, and the 
nearby asteroids Pallas and Egeria when CIT6 was unavail- 
able. The standards have well-tabulated 350 /im flux densi- 
ties, which are available from the CSO/SHARC-II calibra- 
tion web pag^E The final calibration is expected to be better 

1 http : //www. submm. caltech. edu/^sharc 



than 15 per cent, with systematic effects being negligible at 
that level (CSO staff, private communication). 



3 DATA REDUCTION 

3.1 Map-Making and Source Extraction 

We make maps from the raw data and then extract 350 /im 
fluxes from the maps. 

The map- making data reduction package is CRUSH (Com- 
prehensive Reduct ion Utili t y for SHARC-II) , a Java-based 
tool developed by iKovacs! (2006). The software iterates a 
least-squares algorithm to solve for celestial emission along 
with instrumental and atmospheric contributions to the to- 
tal power signal. CRUSH accesses the 7350 Mm polynomial fits 
to obtain a low noise 350 /im sky opacity-based estimate of 
the atmospheric signal. Using the 'deep' utility, as is recom- 
mended for sources fainter than 100 mjy, the maps of each 
field are coadded on a grid of 1.62 arcsec square pixels. The 
outer four rows of pixels are removed automatically by CRUSH 
since they are not sampled sufficiently well to converge to 
useful measurements. The data are fitted with a single Gaus- 
sian beam profile with a FWHM of 9 arcsec. Thumbnails of 
the reduced maps centred on the SHADES 850 /im positions 
are shown in Fig.[T] The bulk of the structure in these images 
is detector noise. 



3.2 Tests of the Data 

Peaks are identified in maps of the signal-to-noise ratio 
(S/N) using an algorithm which only accepts sources sep- 
arated by at least 3<Tbeam = 11.5 arcsec. The maps are also 
inverted and negative sources are extracted in the same way 
as for the positive map. The numbers of peaks found above 
a given threshold are listed in Table as are the numbers 
of peaks expected in a Gaussian field. The total area of the 
maps searched for peaks is 44.5 arcmin 2 . 

Notice that there are more negative peaks detected at 
any given noise level than are expected under the assump- 
tion that the variance is what one calculates from the time 
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Figure 1. 30 X 30 arcsec 2 SHARC-II 350 fim thumbnail images of each SHADES source, centred on the 850 /an positions and using 
1.62 arcsec pixels. Each thumbnail is extracted from a map that has been convolved with the FWHM=9 arcsec SHARC-II point spread 
function. The images are all displayed on the same flux density scale for comparison. LOCK850.6, 63 and 77 lie on the map edges. See 
Tablc[3]for the corresponding flux densities, noise estimates, S/N values, and counterpart associations. 



stream. However, if the maps have Gaussian distributed 
noise, but with a variance 7 per cent larger than we infer 
from the time series the Negative and Gaussian columns in 
Table [5] would match; the number of Gaussian peaks ex- 
pected at S/N = 0.93 x (2,2.5,3) is (67.0, 25.0, 6.3). We 



find this level of agreement between variance in the time 
series and variance in the maps encouraging, and use the 
number of negative peaks to set the statistical significance 
of any detections of positive flux in the maps. 
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Table 2. Numbers of positive and negative peaks in S/N maps. 





2<T 


2.5 o- 


3<T 


4 CT 


Positive Peaks 


72 


32 


9 


1 


Negative Peaks 


66 


28 


7 





P N (10") 


12.9% 


5.4% 


1.4% 




Gaussian Model 


51 


16 


3.4 


0.6 


P N (10") 


10.0% 


3.1% 


0.7% 





3.3 350 /im flux densities of SHARC-II-detected 
SMGs 

If a statistically significant peak is found in the SHARC- 
II map within a small search radius (discussed below) of 
a SHADES position, we infer that this is the counter- 
part source and ca n use its flux, after deboosting as per 
ICoppin et all (|2005h . to constrain the source SED. 

There are several different approaches to estimating the 
350 /im flux density associated with SHADES sources if a 
SHARC-II source is not detected, and they are biased in 
different ways. (1) Measuring the flux density of each ob- 
ject at the 850 /im position will be biased low on average 
because of the uncertainty in the true source location due 
to the large beam sizes and modest S/N in the SHADES 
maps. (2) Measuring the flux density of the brightest pixel 
within given radius of the 8 50 /xm position wil l be biased 
sightly high on averag^3 (e.g. ICoppin et al.|[2005l ; see discus- 
sion below). (3) Measuring the 350 /jm flux density of each 
object at the precise radio position, if one exists, is less bi- 
ased on average than methods 1 and 2 because of the small 
positional uncertainty. 

Not all of the 850 /jm sources mapped at 350 /im have 
radio counterparts. We therefore choose to measure the flux 
densities using method 2 and to correct for flux boosting 
effects, since this measurement can be performed for all of 
our target SMGs in a uniform way. Detections, and measure- 
ments from methods 1, 2 and 3 are given in Table[3] without 
deboosting, for completeness and inter-comparison. The de- 
boosted fluxes are listed in Table [4] In fact we deboost the 
detected sources and the other sources in the same way. 

The first step in either source detection or method 2 is 
to choose a search radius for companion sources. Too large 
a radius increases the false detection rate and the flux boost 
factor. Too small a radius causes sources to be missed. The 
rows labeled Pn(10") in Table[2]show the percentage chance 
that an arbitrary 10 arcsec circle contains a peak. Using the 
negative peaks as a reliable measurement of the noise level in 
the maps, we infer that a 2.5 a peak found within 10 arcsec of 
a SHADES location is 5 per cent likely to be a false positive 
association, i.e. a 2.5 a peak found within a 10 arcsec radius 
is a 95 per cent confidence detection. 

We use Monte Carlo techniques on the actual data to 
test the number of false identifications made. An area with 
a 10 arcsec radius around each 850 /im source is masked out 
so that any real counterpart will not contaminate the test. 
We select a random position on each map and search for a 

2 llKovacs et alj l200rj) and lLaurent et"al] |2006| adopt this ap- 
proach for measuring fluxes for detections and non-detections 
alike. 



peak in the map within the given search radius above the 
designated S/N threshold. This is repeated 10,000 times over 
the SHARC-II maps, and we find that using a search radius 
of ~ 10 arcsec and S/N> 2.5 finds a SHARC-II peak at 
random in 5 per cent of the trials, confirming the conclusion 
above from Tabled We therefore adopt counterpart search 
criteria of 2.5 a and 10 arcsec. 

We can estimate the fraction of sources for which a 
10 arcsec search misses the true companion due to po- 
sitional uncertainty i n both the SHADES and SHARC- 
II data. Ilvison et all (|2007l ) find a one-dimensional po- 
sitional uncertainty of 3.2 arcsec for SHADES SMGs 
by comparing the 850 /jm determined positions (beam 
FWHM=14.8 arcsec) with more precise radio positions (syn- 
thesised beam FWHM=1. 3 arcsec). This is consistent with 
th e theoretical expec tation (see equation 2 of Appendix B 
in Ilvison et all |2007f > of a ~ 0.6 (S/N)" 1 FWHM ~ 3 for 
FWHM=14.8 arcsec and S/N ~ 3. The positional uncer- 
tainty of SHARC-II observations stems from the telescope 
pointing uncertainty, typically 2-3 arcsec. Adding these un- 
certainties in quadrature and converting one-dimensional 
uncertainties to two dimensions and then a radial distance, 
we expect the true 350 /im counterpart to lie within a 
10 arcsec circle 93 per cent of the time. Using a larger search 
radius would increase the level of flux boosting without in- 
cluding more sources. Peaks in the SHARC-II map will be 
displaced from the true counterpart location by a distance 
that scales with the SHARC-II beam size (9 arcsec) and the 
S/N (> 2.5a), as for SHADES. Adding this term, we find 
that 90% of apparent peaks associated with a counterpart 
will lie within our chosen 10 arcsec circle. 

Proper correction for the flux boosting that results from 
picking off peaks in low S/N data requires knowl edge of the 
source count distribution (e.g. lCoppin et al]|2005l ). However, 
prior information about the 350 /mi source counts is not 
yet sufficiently well-constrained. We have estimated the flux 
boosting effect on this sample of low S/N 350 /im-detected 
and non-detected SMGs (< 4 a) by calculating the error- 
weighted mean excess flux of method 2 compared to method 
3 to be a factor of ~ 1.5. We therefore divide each 350 )im 
flux measurement, from method 2 (or from modest S/N de- 
tections), by 1.5 and present the deboosted fluxes in Table|4] 
for use in fittin g SEDs. For the high est S/N 350 /im-detected 
SMGs (> 4.5 q; lKovacs et alj |2006). flux boosting effects are 
negligible and so we do not correct these fluxes. 

3.4 Discussion of other 350 /im photometry of 
850 /im SHADES sources 

The 350 ytm flux densities have been constrained for 25 per 
cent of the SHADES 850 /jm sources, with the current work 
providing 21/31 of the SHADES 350 /im flux density con- 
straints. Complementary SHARC-II observations of SMGs 
including seve r al SH ADE S sources have been performed by 
lLaurent et "ail (|2006T ) and lKovacs et all (|2006t ). Since we in- 
clude their photometry in the SED fits of seven SHADES 
sources, it is worth describing how fluxes are measured by 

thes e groups. 

iKovacs et al.l (|2006T ) conducted follow-up observa- 
tions of bright (>5mJy) radio-identified SCUBA sources 
with optical spectroscopic redshifts, including LOCK850.3, 
LOCK850.14, LOCK850.18, and LOCK850.30. They de- 
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tected 12/15 of their sample with a mean S/N ~ 4.5 for the 

det ections. 

lLaurent et alJ (|2006l ) performed follow-up observa- 
tions of 17 Bolocam 1.1 mm-selected source candidates 
in the LH with SHARC-II. Of the 17, ten are de- 
tected, including LOCK850.1, LOCK850.2, LOCK850.3, 
LOCK850.12, LOCK8 50.14, LOCK850.41 LOCK850.27, 
and LOCK850.7fB (see lLaurent et alJliooBI for counterpart 

ide ntification j3) ■ 

iKovacs et al.l (|2006T l and lLaurent et al] ||2006T ) quote as 
a flux the peak flux density in an specified search radius 
around their target source position, which is taken from ei- 
ther radio or mm data. They do not deboost these fluxes, 
but their observations are at lower noise levels so this step 
is less important. 

Three SHADES sources that we observe in this pro- 
gramme, LOCK850.1, LOCK85 0.3 and LOCK850.41 have 
also been be e n ob served by IKovacs et all (|2006l ) and 
lLaurent et ail (2006). Because their observations are at 
lower noise we use their fluxes in Tableland in subsequent 
analysis. 



" 




Figure 2. A multi- wavelength view of LOCK350.1, a 
3.9 cr 350 /im blank-field source detected at RA=10 h 52 m 43?2, 
Dec.=57°23'9 / ' (J2000). The 30arcsec X 30arcsec images are 
centred on the position of LOCK350.1. Left: High-resolution 
(1.3 arcsec) 1.4 GHz c ontours (-3, 3, 4.. .10, 20.. .100 X a) from 
Bigg s fc Ivisonl l|2006h overlaid on greyscale archival i?-band opti- 
cal data. Right: Low-resolutio n (4.2arcsec) 1.4 GHz c ontours (—3, 
3, 4.. .10, 20. ..100 X cr) fr om iBiggs fc Ivisonl 11200611 overlaid on 
greyscale 24 /im data from llvison et al. I J2007T) . 



3.5 Formally detected 350 /mi counterparts 

In the analysis section of this paper we use our best estimates 
of the 350 /im flux of each source as photometric constraints. 
However, for comparison to blank field surveys it is useful to 
consider how many of the sources are seen at a high enough 
S/N that they would be counted as detections. 

We have identified seven 350 /im counterparts of 850 /mi 
sources, as presented in Table [3] The positional offsets of 
the detections from the SHADES positions are consistent 
with the positional distribution discussed above; 5/7 (71 per 
cent) of the detections lie within 5 arcsec. No trend is seen 
of decreasing positional offset with increasing S/N ratio of 
the seven associations, but this is not surprising given the 
low number statistics and the small dynamic range in S/N. 

It is striking that the number of detected counterparts 
to SHADES sources exactly matches the excess of positive 
peaks over noise (or negative peaks) in the full mapped area, 
as listed in Table [2] 

As we did with the number of peaks, we can use the cor- 
responding composite negative catalogue to test if real as- 
sociations are being found or just noisy peaks in the 350 /mi 
maps. Counterparts to 850 /mi sources are searched for in an 
inverted map in the same way, using a S/N threshold of 2.5 
and search radius of 10 arcsec and no 'negative' counterparts 
are identified. This result reassures us that the 350/850 /im 
associations are likely to be real. 



3.6 350 /mi flux density constraints of SMGs not 
individually detected by SHARC-II and flux 
boosting 

Any of the following mechanisms might cause an SMC not 
to be detected at 350 /im: (1) the 850 /im source could be 
intrinsically less luminous, colder and/or have a different 
SED shape than the strongly detected sources; (2) the SMG 
could lie at a redshift in excess of about 3, this being more 
likely if there is also no radio counterpart, since the existing 
radio data are sufficiently deep to obtain a l arge fraction 
of co unterparts for SMGs only for z < 3 (see llvison et al.l 
120051 ) . and above z ~ 3, the 350 /im band is sampling the 
Wien side of the SED and suffers from cosmological dim- 
ming without the benefit of the negative K-correction, or 
(3) The SMG could be spurious, although the false positive 
rate in the SHADES catalogue is believed to be very low. 
The catalogue was carefully constructed from a comparison 
of multipl e data reductions in order to minimise false detec- 
tions ("see lCoppin et al.ll2006T l. 

We rely on the photometric fluxes even for the 16 
sources which are not detected at high significance. There- 
fore we have performed a stacking analysis to determine if a 
significantly positive 350 /im flux density is associated with 
those sources. 

350 /im flux densities and errors are measured on the 
SHARC-II maps at the radio position for the 11 of these 16 
sources which have a radio counterpart and we find the mean 
flux is 16.7 ± 4.8 mjy. This is our least biased flux estimator, 
and the agreement with our listed deboosted fluxes is very 
encouraging. The weighted mean deboosted flux for all 16 
of these low significance sources is 16.6 ± 3.8 mjy. 



3 LOCK850.3 and LOCK8 50.14 are detected b y IKovacs et all 
II2006T ) but also discussed by lLaurent et aL (2006). 

4 LOCK850.27 and LOCK850.76 were not previously associated 
with Bolocam sources, since they are new SHADES detections. 
LOCK850.27 and LOCK850.76 lie 4.8 and 10.6 arcsec away from 
two Bolocam sources, and so two new SCUBA/Bolocam associa- 
tions are claimed here. 



3.7 Serendipitous 350 /mi blank-field SMGs 

We identify one serendipitous 3.9 a 350 /im source in our 
map of LOCK26/32: LOCK350.1. LOCK350.1 is not as- 
sociated with any 850 /im sources, since it is > 15 arcsec 
away from any of the SHADES 850 /im positions in the 
map. Based on Gaussian noise statistics and the number 
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Table 3. SHARC-II 350 /an flux densities of 850 /an sources from SHADES. Coordinates are only given for > 2.5 a SHARC-II detections 
< 10 arcsec away from the SHADES 850 /im position. Each offset has an uncertainty of about ±2 arcsec, reflecting the uncertainty in 
the SHARC-II map astrometry due to the RMS pointing uncertainty. The flux densities in the 'brightest pixel' column represent the 
highest S/N pixel in the beam-convolved image within 10 arcsec of the SHADES catalogue position - these fluxes are used in the SED 
fitting once a correction has been made for flux boosting (see Table|4]for the corrected fluxes). Non-detections are given in parentheses; 
measurements were also made on the SHARC-II maps at the SHADES 850 fim and radio positions for comparison. See the text for a 
discussion of the biases of each of these measurements. The observations for LOCK850.15 are insufficiently deep (incomplete) and its 
map is thus left out of the general analysis (see text). 



SHADES ID 


350 /im position 




Offset 




^350 




Comments 




RA 


Dec. 






brightest pixel (< 10 ) 


at 850 |im posn. 


| at radio posn. 






(J2000) 


(J2000) 




(arcsec) 


(mJy) 


(mJy) 


(mJy) 




LOCK850.01 








(~ 4.3) 


(23.0 ± 21.2) 


12.7 i 22.2 


22.6 i 21.4 


see Laurent et al. (2006*1 


LOCK850.03 


1 fl 11 c i9 m 38 s 7(1 


57° 24' 37'' 


4 


3.7 


67.1 ± 18.5 (3.6 a) 


59.7 ± 18.0 


65.0 ± 18.2 


see Kovacs et al. (2000, 




















LOCK850.04 


10 h 52 m 04?42 


57° 26'54(' 


3 


5.0 


37.3 ±9.1 (4.1 a) 


21.0 ±9.2 


18.9 ± 9.3 


















no rj _|_ n "1 

oo.Z ±9.1 




ljL>Oi\.8oU.UD 


- 


- 




(~ 10) 


(D 1 .U ± O t .0) 


1.8 ± 43. 1 


2.6. D ± J9.1 


cage or briArtO-11 map 










(~ 10) 


(3D. 2 ± 12.0) 


1 1 11 J_ 11 fl 

14. y ± 11.9 


1 n a _i_ too 

iy.4 ± lz.z 




LOCK850.15 








(~ 4.2) 


(9.8 ± 55.7) 


-32.7 ± 53.7 


-8.3 ± 52.5 


incomplete observations 
















6.9 ± 54.9 




LOCK850.16 








(~ 4.9) 


(38.6 ± 15.8) 


26.6 ± 15.5 


33.3 ± 15.6 




LOCK850.21 








(-6.5) 


(25.1 ± 13.2) 


14.3 ± 12.9 






LOCK850.22 








(~ 10) 


(13.1 ± 15.6) 


-1.3 ± 15.6 






LOCK850.26 








(~ 4.5) 


(18.3 ± 8.8) 


6.9 ± 8.5 


13.4 ± 8.6 




LOCK850.28 


10 h 52 m 57=86 


57° 30' 59" 


7 


10.0 


34.9 ± 11.7 (3. Oct) 


8.3 ± 11.8 


(29.6 ± 11.6) 




LOCK850.33 


10 h 51 m 55?82 


57°23'll" 


3 


1.4 


24.7 ± 8.4 (2.9 ct) 


21.5 ± 8.4 


20.4 ± 8.4 




LOCK850.41 








(~ 7.5) 


(42.5 ± 25.4) 


15.9 ± 24.4 


12.4 ± 24.2 


see Laurent et al. (2006") 
















34.8 ± 24.8 




LOCK850.47 








(~ 9.4) 


(24.4 ± 20.9) 


-8.3 ± 21.2 






LOCK850.48 








(-6.5) 


(24.3 ± 13.5) 


9.9 ± 12.3 






LOCK850.63 








(~ 10) 


(52.3 ± 29.1) 


9.9 ± 31.4 


15.6 ± 29.1 




LOCK850.64 








(~ 10) 


(17.4 ± 12.4) 


-5.1 ± 12.8 






LOCK850.77 


10 h 51 m 56?22 


57°22'09" 


8 


6.3 


62.1 ± 24.6 (2.5 <r) 


9.9 ± 18.7 


-0.2 ± 17.9 


edge of SHARC-II map 


SXDF850.1 








(~ 6.2) 


(24.5 ± 13.9) 


15.9 ± 14.5 


15.9 ± 14.6 




SXDF850.3 


02 h 17 m 41?95 


-04°56'26 


'3 


3.4 


39.3 ± 14.3 (2.7<r) 


23.1 ± 14.1 


23.1 ± 14.1 




SXDF850.8 








(~ 6.3) 


(24.5 ± 14.1) 


10.5 ± 14.1 


18.1 ± 14.5 




SXDF850.il 


02 h 17 m 24?81 


-04°59'37 


'2 


4.6 


46.6 ± 15.0 (3.1 <r) 


26.6 ± 14.2 


22.2 ± 14.1 




SXDF850.17 








(~ 10) 


(18.6 ± 21.0) 


-13.0 ± 17.3 






SXDF850.119 








(~ 7.6) 


(23.8 ± 14.9) 


13.0 ± 15.3 


0.8 ± 15.4 





of independent beam sizes in the survey area with noise 
less than 20mJj|3, about 0.1 false positives are expected on 
average at a S/N of 3.9; LOCK350.1 is therefore likely a 
real blank-field 350 /im source. The position and flux den- 
sity of LOCK350.1 are RA=10 h 52 m 43 , :2, Dec -57°23'09" 
(J2000) and 32.8 ± 8.3 mJy. This position corresponds to a 
region of positive flux density in the 850 fj,m map and a 2.3 a 
peak is located ~ 4 arcsec from the position of LOCK350.1. 
There also appears to be a 4 a 1.4 GHz source and a faint 
24 counterpart a t these coordinates jlvison et a.1.1 120(371 : 
iBiggs fc Ivisonl 120061 : see Fig. H}. This source is therefore 
probably the thir d se cure 350 fim blan k-field detection (see 
iKhan et al1l2005l and lKhan et al.ll2007l for the first two). 



4 ANALYSIS AND RESULTS 

4.1 Constraints on the 350 /im source counts 

We now attempt to crudely estimate the 350 /im source 
counts at our typical observed flux density, even though 
the SHARC-II data are not blank-field observations. We 

5 A noise cut of 20 mjy was chosen arbitrarily, since higher 
S/ N detections in noisie r regions are more likely to be spurious; 



place an upper limit on the source counts using the num- 
ber of detections acquired in the observed area; one would 
expect to do worse than this in a blank-field survey, since 
here known SMGs were observed. Given a total observed 
area of 48.3 arcmin 2 and seven detections above a 350 /im 
flux density of ~ 25mJy, we estimate N(>S) < 500deg~ . 
We estimate a lower limit on the source counts by apply- 
ing the SHARC-II detection rate of 850 /im sources to the 
whole SHADES area; one would expect to do at least as 
well in a blank-field search. We estimate a lower limit of 
JV(> S) > 200 deg" 2 , given a SHARC-II detection success 
rate of 7/24 and a total of 120 SHADES sources found in 
~ 720 arcmin 2 . Note that we have neglected the possibil- 
ity of enhanced density from clustering. The limits derived 
here a re consistent with the predictions from lLagache et al.1 
J2004h and scenario E from the semi-analytic models of 
iGuiderdoni et all (|l998l ) (see Fig. |3). Although crude, this 
is the best available es timate of th e 350 source counts 
at these flux densities. IKhan et al.l l|2007l ) also report the 
350 /im source counts above 13 mJy from a deep survey with 
SHARC-II. Large surveys planned with the Balloon-borne 



Large Aperture Submm Telesc ope (BLAST; Devlin et al 
2004), Herschel (iPilbrattJ I2003T). SCUBA- 2 |Holland et al 



2006), and Planck l|Tauberl 12004) will be able to provide 



further constraints on the source counts at shorter submm 



e.g. ICoppin et all J2005h . 
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SHADES ID S350 S S 5o S450 Si.4GHz S24 Si.i mm Si. 2mm phot z spec z 

(mjy) (mjy) (/*Jy) (^Jy) (mjy) (mjy) 



LOCK850.04 


24.9 


± 9.1 


10.65 


(±11) 

V- 1 - 1 . 8 / 


< 


134 




32.0 ± 5.1 
73.0 ± 5.0 


201 
179 


± 73 
± 68 


LOCK850.06 


38.0 


± 37.6 


6.85 


(±H) 


< 


77 


15 


± 4.8 (tent.) 


75.1 


± 12.7 


LOCK850.10 


24.1 


± 12.0 


9.15 


(±11) 


< 365 




25.5 ± 6.3 






LOCK850.15 


6.5 


± 55.7 


13.25 


(±ln) 


< 


14!) 




43.9 ± 7.8 


353 


± 20 
















61.5 ± 7.6 






LOCK850.16 


25.7 


± 15.8 


5.85 


(±M) 


< 


67 




106.0 ± 6.0 


314 


± 24 


LOCK850.21 


16.7 


± 13.2 


4.15 


(±l:f) 


< 


70 




5 a < 30 


97.9 


± 14.1 


LOCK850.22 


8.7 


± 15.6 


7.55 




< 


76 




5 a < 30 


402 


± 21 


LOCK850.26 


12.2 


± 8.8 


5.85 


(±1:1) 


< 


48 




31.4 ± 5.2 


195 


± 16 


LOCK850.28 


23.3 


± 11.7 


6.45 


(±1) 


< 


56 










LOCK850.33 


16.5 


± 8.4 


3.85 


(±i:?) 


< 


49 




51.0 ± 4.3 






LOCK850.47 


16.3 


± 20.9 


3.55 


(± 2 :I) 


< 


21 










LOCK850.48 


16.2 


± 13.5 


5.45 


(±11) 


< 


79 






203 


± 17 


LOCK850.63 


34.9 


± 29.1 


3.65 


(±1:1) 


< 


50 




22.6 ± 4.8 


236 


± 17 


LOCK850.64 


11.6 


± 12.4 


5.85 


(±1:1) 
(4:1) 


< 


95 










LOCK850.77 


41.4 


± 24.6 


3.25 


< 


39 


15 


5 ± 4.4 (tent.) 


51.7 


± 13.1 


SXDF850.1 


16.3 


± 13.9 


10.45 


(±1:1) 


< 


65 




54.3 ± 9.7 






SXDF850.3 


26.2 


± 14.3 


8.75 


(±ii) 


< 


81 




77.2 ± 9.3 


110 ID 


SXDF850.8 


16.3 


± 14.1 


6.05 


(±1:1) 


< 


45 




52.0 ± 9.5 






SXDF850.il 


31.1 


± 15.0 


4.55 


(± 1 2-b 


< 


79 




50.8 ± 10.0 


195 


± 47 


SXDF850.17 


12.4 


± 21.0 


7.65 


(±li) 
(±lh 


< 


71 










SXDF850.119 


15.9 


± 14.9 


4.55 


< 


70 


38 


± 9.7 (tent.) 


784 


± 47 



(0.526 or 1.482) 2 



1.8 ± 0.5 
1.6 ± 0.4 



1.6 ± 0.4 



1.7 ± 0.4 



3.1 (±£|) 

2 4 (±8:t) 

1.9 (±g;«) 

> 1.0 

> 2.0 
3.6 (±»;i) 

> 2.0 



> 1.5 

2.4 (±°;f) 

2.6 (±8:1) 

> 1.5 
2.6 (±°i) 
2.6 (±lp 
2.1 (±°;f) 
2.6 (±i;f) 

2.4 (±g:l) 

> 2.0 

2.1 (±g : }) 



2.664 2 , 2.686 1 , 3.699 3 



LOCK850 


01 


24.1 


± 5.5 


8.85 (±i:g) 


< 


47 


78.9 


± 


4.7 


217 


± 


10 


4 


4 


± 


1 


3 


3.6 


± 0.5 


2 


4 




I) 


2.148 1 ' 2 


LOCK850 


02 


25.3 


± 10.3 


13.45 (±1;}) 


< 


123 


40.7 
52.4 


± 
± 


5.6 
5.6 


545 


± 


31 


6 


8 


± 


1 


4 


5.7 


± 1.0 


2 
2 


9 
9 




?) 
I) 




LOCK850 


03 


40.5 


± 6.5 


10.95 (±J:|) 


< 


34 


35.0 


± 


5.2 


183 


± 


33 


4 


8 


± 


1 


3 


4.6 


± 0.4 


2 


6 


!±8 


i) 


(2.94 orS.036 1 ) 2 
















25.8 


± 


4.9 


175 


± 


23 
























LOCK850 


12 


29.3 


± 16.0 


6.15 (±J p 


< 


35 


44.3 


± 


5.1 


203 


± 


19 


4 


1 


± 


1 


3 


2.6 


± 0.4 


2 


(i 




I) 


2.142 1 ' 2 


LOCK850 


14 


41.0 


± 6.8 


7.25 (±i;|) 


< 


96 














5 


1 


± 


1 


3 


3.4 


± 0.6 


2 


(i 


(± 


?> 


2.611 1 ' 2 ' 3 


LOCK850 


18 


7.5 


± 6.7 


6.05 (±l'f) 


< 


84 


29.4 


± 


4.4 








5 


1 


± 


1 


3 


3.4 


± 0.6 


3 


1 


(± 


T> 


1.956 1 ' 2 


LOCK850 


27 


3.4 


± 5.1 


5.05 (±M) 


< 


32 














5 


2 


± 


1 


4 


3.2 


± 0.7 


4 


6 


(±0 


1) 




LOCK850 


30 


38.0 


± 7.2 


4.75 (±};|) 


< 


86 


245 


± 


13 


233 


± 


19 












(0.4 


± 0.8) 


2 


1 


(±:: 


I) 


2.689 1 ' 3 


LOCK850 


41 


10.3 


± 5.5 


3.85 (±?:§) 


< 


16 


43.6 


± 


4.7 


051 


± 


46 


4 





± 


1 


3 


2.4 


± 0.5 


3 


4 


(±8 


I) 


0.689 1 ' 2 ' 4 














22.1 


± 


4.8 


475 


± 


37 
























LOCK850 


76 


4.4 


± 6.7 


4.75 (±11) 


< 


90 


48.0 


± 


6.0 


592 


± 


20 


4 


4 


± 


1 


4 






4 


6 


(±1 


i) 





Table 4. Multi-wavelength photometry of SHARC-II-observed SHADES sources. The first set are SHADES sources th at use 350 )im 
flux densities measured here, and the second set are 350 fira flux densities of SHADES sources obtained from the literature llKovacs et al.l 



120061 and lLaurent et al.ll2006T). 350 um p hotometry has been corrected for flux boosting (see Section 13. 6) . 850 (deboosted) and 450 ^tm 
photometry are from lCoppin et al.N2006l) , 1 .4 GHz and 24 (i m dat a are from llvison et alT l l2007l ) ('tent.' refers to tentative identifications), 
1.1mm and 1.2 m m photometry are from lLaurent et al] l|2005l ) and Dunlop (private communication; an improved reduction of the 
iGreve et al 1l2004l MAMBO data), respectively. In cases wh ere two robust radio an d / or 24 /an IDs exist for a single SHA DES source, both 
are listed. Photometric redshifts are from tables 3, 4 and 5 of|Aretxaga_et^il and spectroscopic redshifts are from: 1. Chapman et al. 

1120051 ); 2. llvison etHI 1120051 ); 3. IChapman et al] i2003l) ; 4. IChapman et al.l hOCO) . 



wavelengths and over a wider dynamic range of flux densi- 
ties. 



4.2 SEDs of SMGs 

The SEDs of SMGs is dominated by thermal emission from 
cold dust. SHARC-II photometry of the SHADES sources 
can determine the apparent temperature of the dust. In con- 
junction with knowledge of the redshift, this allows inference 
of dust masses, FIR luminosities, and SFRs, allowing one 
to place these objects in context with other populations of 
high-redshift star-forming galaxies and AGN. 

The shape of a luminous dusty galaxy SED is well ap- 
proximated by a modified blackbody spectrum described by 
the dust temperature, Td, and dust emissivity e oc . where 
P lie s in a physically plausible range of 1-2 (Hi ldebrandl 
1983). There is a degeneracy betwee n Td and (3 which can- 
not be disentangled by our data dBlain et al.l |2002| ). so 
we fix (5 = 1. 5 , whi ch is consistent with the finding of 
Dun ne fc Eales ( 200 ll ) for local galaxies and with the val- 
ues obtained for carbonite and silicate grains from labora- 



tory measurements l|Agladze et al.l Il994 ). This leaves two 
free parameters to be fitted: Td and the SED normalisation. 

In practice we work in terms of the apparent dust tem- 
perature, Ta = Td/(l + z) and fit the data for each galaxy to 
Af^f /[exp(h v bs/k Ta) — 1], where several factors of (1 + z) 
have been absorbed into the apparent normalisation, A. 

Assuming that the redshift is known, the restframe dust 
temperature can be recovered and the luminosity of the 
SMG can be determined by integrating the SED. In cases 
where the spectroscopic redshift is ambiguous, the most reli- 
able or probable spectroscopic redshift available from the lit- 
erature has been chosen and is listed in Table [S] For sources 
with only a 90 per cent confidence photometric redshift lower 
limit, the redshift is taken at this limit (see Table |4}. 

We characterise the shape of the SEDs independently 
of the photometric redshifts because we do not use the radio 
data in our fits to the SED, and where photometric redshifts 
are used they have been determined using all available FIR 
and radio photometry (except for the 350 /im data from this 
paper) . 

More complex SED modelling was not attempted 
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10 5 F 



LO 3 



10" 
0.1 




- - Guiderdoni et al. (1998) A 

Guiderdoni et al. (1998) E 

— Lagache et al. (2004) 
Baughetal. (2005) 

A Khan et al. (2007) 

# This work 



1.0 10.0 
Observed 350 u.m Flux Density S [mJy] 



100.0 



Figure 3. Cumulative 350 (im source counts. The s o lid lin e is 
a prediction of the source counts from lLagache et al. us- 
ing a phenomenological model that constrains the IR luminos- 
ity function evolution with redshift. The dot-dash line is the 
pr ediction of the 350 counts from the semi-analytic models 
of iBaugh et al. I 1120051 ) which match the 850 /im number counts. 
The dotted and dashed lines are semi-analytic predictions of the 
source counts re sulting from evolu t ionary scenarios E and A, re- 
spectively, from Guiderdo ni et aTl (jl99ST) ■ where star-formation 
proceeds in 'burst' mode and includes either no ULIRGs (sce- 
nario A) or an increasing fraction of ULIRGs with redshift (sce- 
nario E). The solid circles with arrows are upper and lower lim- 
its on the source counts based on this work (see text) and are 
consis tent with the sourc e count predictions of [Lagache et al.l 
[120041) and scenario E of iGuiderdoni et all lll998ri ~ scenario A 
of IGuiderdoni et al. is inconsistent with the observational 

limits, as expect ed since it contains no ULIRGs. The source count 
constraint from Khan et al. above 13 mJy is also plotted 

(open triangle). 



(e.g., fitting two-te mperature dust components as in 
iDunne fc E alcs 2001), since this would require more free 
parameters than can be constrained using the typically 2- 
3 photometric points which exist for each of our SMGs. 
The Wien side of the spectrum is sometimes also mod- 
ified by a power-law of the form S v <x v~ a to account 
for the increase in optical depth in this part of the spec- 
trum and to provide a better fit t o observational data (see 
iBlain. Barnard fc Chapman! l2003r). We ne glect such elabo- 



rations here (see also Kovacs et alj [200(3 ) . since the Wien 
side of the spectrum is not sampled with our data. 24 /im 
photometry are available for SHADES sources but are com- 
plicated to interpret since this band samples PAH and stel- 
lar e mission features in high-re dshift star-forming galaxies 
(e.g. ISaiina. Lacv fc Scotdf2005l ). 

The results of this fitting process are plotted in Fig. U 
Notice that the 350 /jm data point provides almost all of the 
contraining power on Ta. 

We have measured the apparent temperatures (Ta) of 



an unbiased sample of SMGs from SHADES. We calcu- 
late the corresponding temperature distribution by adding a 
rectangle of area one centred on the source Ta, with a width 
given by the errors in Ta to the histogram (see Fig. [5}. It 
appears that our survey selects SMGs with Ta — (8 ± 3) K. 
We recover the intrinsic dust temperature (Td) distribution 
of SMGs by shifting each SED into the restframe, follow- 
ing a similar procedure to that described for the Ta dis- 
tribution above (see Fig. [SJ). The Td distribution appears 
to contain structure ('bumps' centred around 12 K, 28 K, 
and 35 K) and is broader compared with the Ta distribu- 
tion. To investigate if this distribution is more structured 
and/or wider than we would expect by chance, we have sim- 
ulated the temperature distribution we would expect to see 
if all SMGs have T d = 28 K. The simulated T d distribu- 
tion is constructed by adding an error box of area one cen- 
tred on Ti = 28 K + rand x a, where 'rand' is a normally- 
distributed random number, for each source. Although the 
shape and width of the simulated distribution changes no- 
tably depending on the random realisation, we never re- 
produce the 'bumps' at Ti ~ 12 or 35 K and we always 
predict more SMGs at a few K than we see. Therefore we 
conclude that the underlying distribution is just a bit wider 
than we expect by chance and is consistent with most SMGs 
at 28 K with a few hotter ones. Now we compare the dust 
temperatures of our photometric and spectroscopic redshift 
sub-samples and to previous estimates for SMGs. 

We derive a mean T d = (31 ± 18) K (median T d = 31 K) 
for our sample of SMGs with spectroscopic redshifts. For 
the sample of 21 350 ^im-observed SMGs with only photo- 
metric redshifts, the mean is Td = (35 ± 17) K (median 
Ti = 29 K). A two-sided KS test reveals that the temper- 
ature distributions of the photometric redshift and spec- 
troscopic redshift subsets are consistent with being drawn 
from the same distribution. Therefore we quote a mean 
T d = (34 ± 17) K (median T d = 29 K) for our entire sample 
of SMGs. Our mean Td is consistent w ith previous estimates 
for SMGs from lChapman et all l|2005h (T d = 3 6 K for radio- 
identi fied SMGs spectroscopi c redshifts), f rom | Kovacs et ail 
|2006l ) (T d = 35 K), and from lPope et""ai] (|2006l ) (Td ~ 30 K 
for SMGs in the Great Observatories Origins Deep Survey- 
North, GOODS-N , region). In contrast, we find that the 
iFarrah et aT]l2003l ULIRG sample has a median T d = 42 K 
when their photomet ry are fit with our greybody (see also 
IClements et al. 1 120081 ). Thus SMGs appear to have lower 
temperatures on average than local ULIRGs, although bear 
in mind that local ULIRGs are selected at 60 /im which in- 
herently introduces a bias to selecting hotter objects than 
submm observations would typically find. 

The SED is integrated to estimate the bolometric lu- 
minosity of the SMG. Assuming that the FIR luminosity 
is predominantly powered by star-formation (i.e. negligible 
contribution from an AGN), coming f rom a starburst of less 
than 100 Myr, with a ISalpeterl l|l955l ) ini tial mass function , 
then the SFR can be calculated following iKennicuttl l|l998l ): 



SFR(Mq yr _1 ) = 1.7 x 10~ 10 L F m(L c . 



(1) 



The median FIR luminosity for our sample is ~ 2 x 
1O 12 L (the mean is ~ 8 x 10 12 Lq , with a large scatter) , 
which is consis tent with what has been found previously 
for SMGs (e.g. iKovacs et aLI l2006l l . The median SFR for 
the 350 /tin-observed SHADES sources is ~ 400 M yr _1 , as 
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Figure 4. Best-fitting modified blackbody SEDs for SHADES 850 /xm-selected SMGs with 350 /im photometry with fixed to 1.5, 
in the same order as Table [5] The first ten SEDs use the available spectroscopic redshifts, and the rest use only photometric redshift 
constraints. The dashed curves show the 1 a range in the best-fit SED. The solid circles indicate photometry used to fit the SED, while 
the open circles show the SCUBA 450 ^m upper limits which are not used explicitly in the fits (though they are never violated). The 
spectroscopic or photometric redshifts indicated in each panel have been used to recover the given restframe Td . Notice that the 350 fim 
point provides almost all of the constraining power on the dust temperature. 



calculated from Equation [T] This is consistent with previous 
work, indicating that SMGs are significant contributors in 
the global star-formation of the high-redshift Universe (see 
e.g.lLillv et al.lll999l;|Blain et al.ll2002l ; IChapman et al.ll2005l ; 
lAretxaga et al.ll2007f ).~ 

We can also constrain the cold dust mass, Md, in these 
obscured star-forming galaxies. The flux density of a galaxy 



at a n observed frequency, ^ Q b s , is given by the usual relation 
(e.g. iHughes. Dunlop fe Rawlingslll997h : 

S„ oba = B v , (T)*v M d (1 + z)/Dt, (2) 

where Dl is the cosmological luminosity distance and B v i 
is the Planck function evaluated at the emitted frequency, 
v' — Vobs (1 + z). The quantity ka is the wavelength- 
dependent mass-absorption coefficient (or 'effective area' for 



© 0000 RAS, MNRAS 000, 000-000 



SHARC-II follow-up of SHADES sources 1 1 



Table 5. Derived dust temperatures (Td), FIR luminosities (Lfir), illuminated dust masses (Mj) and SFRs for the 350 /mi-observed 
SHADES SMGs. The emissivity, j3, is fixed at 1.5 and the best spectroscopic or photometric redshifts are used in fitting to the 350 
and 850 fim photometry, plus 1.1 and 1.2 mm data if available. Errors in the fit parameters (i.e. the position of the peak and the SED 
normalisation) are derived from the 68 per cent \ 2 confidence interval around the best-fitting values and have been carried through to 
all the derived quantities. Where appropriate, errors in the photometric redshifts have also been folded into the uncertainties. 



SHADES ID 


z 


T d 


-^FIR 


M d 


SFR 






(K) 


(X1O 12 L ) 


(X10 9 M ) 


(Mgyr" 1 ) 


With a spectroscopic rcdshift 












LOCK850.01 


2.148 


27^ 


, ,+ 1.0 

■^-o.s 




, qn +170 
Jau -130 


LOCK850.03 


3.036 


39± 4 4 


9 - 6+ 2's 


o.a±g : » 


moot™ 


LOCK850.04 


0.526 


13+ 2 2 




, 7 +0.9 
^•'-0.6 


ia±S 


LOCK850.12 


2.142 


33+ 14 


3 ,+8.2 


7+ 1 ' 4 
u ''-0.3 


5S0 +1400 
aav ~440 


LOCK850.14 


2.611 


39+1; 


s. +_tt 


5+°- 2 
u -°-o.i 


1400+ 790 


LOCK850.16 


1.147 


25+? 




9 +0S 


140+ 280 


LOCK850.18 


1.956 


18+^ 


5+ - 6 


, 7 +17.0 
^•'-1.3 


„ +ioo 

a -70 


LOCK850.30 


2.689 


76 + 103 


4 „ q+932.0 
* uu -26.0 


1+ ' 1 


+ 101)1)1)1) 

D»UU_ 4500 


LOCK850.33 


2.664 


31 + 8 
' 3i -13 


1 9 + 21 
la -1.5 


5+ 2 ' 1 
u.o_ 2 


33„+370 


LOCK850.41 


0.689 


12+ 4 
1 -6 


04+ 05 


, 7 +12.0 
1 -0.8 


7+9 
'-6 



Kovacs 


ct 


al. 


(2006) 


: Kovacs 


ct 


al. 


(2006) 


Kovacs 


ct 


al. 


(2006) 


Kovacs 


ct 


al. 


(2006) 


Kovacs 


ct 


al. 


(2006) 


Kovacs 


ct 


al. 


(2006) 


: Kovacs 


ct 


al. 


(2006) 



With a photometric rcdshift 

LOCK850.02 

LOCK850.06 

LOCK850.10 

LOCK850.15 

LOCK850.21 

LOCK850.22 

LOCK850.26 

LOCK850.27 

LOCK850.28 

LOCK850.47 

LOCK850.48 

LOCK850.63 

LOCK850.64 

LOCK850.76 

LOCK850.77 



2 9+°' 7 
^•-0.1 

3 6 +1 '° 
, ,+0.9 
•'■ 1 -0.3 

2 4+°' 4 

> 1.0 

> 2.0 

3 6+ - 1 

4 6+ 1 ' 4 

> 2.0 

> 1.5 
2.4j 
2.6 4 



4 

2.6 



+ 0.5 
0.1 
+ 0.4 
0.4 
> 1.5 
+ 1.4 
-1.1 
+ 0.8 
0.1 



+ 12 
-8 
K+51 



30 

55^ 
o 5 +17 
JO -10 
33; 



, + 17 
-29 



> 21 



+8 



-7 



+ 11 
- 15 
7+22 



35 
27_ 

> 29 

> 27 
35 
65 

> 22 

29+ 21 " 
^ a -19 
„,+ 140 
yJ -75 



7 

+ 19 
-19 

+ 20 

-12 

+ 71 
59 
+8 



4 8+ 8 - 3 
^■°-2.9 

19.0+f 9 ° ° 

™±ii 

> o.3t° i 

7+I8.O 

, O + 30.0 
a -°-3.1 

1 9+ 1 ' 2 

> a-o±i'J 

- n „ + 41.0 

> 0.8__ g 

2 6+ 9 ' 4 
^ D -2.0 

,, „ + 2600. 
zz - u -22.0 

> 0.6±g;» 

s, n + 13000 n 
°" i - u -82.0 



, q+2 
!- a -1.0 

3+ 2 ' 5 
U.3_ 2 

9+ 1 ' 2 
ua -0.3 

1 5+ 9 '° 

> o.8±S;| 

> 2.9tf 8 

5+ 1 ' 7 

1 2+ 9 ' 1 
1 -^-1.0 

> 9+°- 9 
5+37.0 

5+ ' 8 
, +850.0 
-0.08 
, + 2.8 
-0.6 
g + 100.0 
ua -0.4 

o.06±°; 4 



0.1: 
> 1.3^ 



„„ n +1430 
=30-490 
,, nn +100000 
MUU -3200 
g 20 +780 
ai: -610 
ego +6400 
Oau -630 



+ 100 
-40 
-3000 
-110 
+ 5200 

)40 
+ 1200 
210 
340+220 

+ 7100 
140 
+ 1600 
-340 
100000 
3800 



> 60 
> 120 
650' 
330 



> 140 

440 
3800 



> 90±^° 

410l|™ 

14000+ 2390000 
14000 



X < 



< 1 



1. unphysical, misID? 



SXDF850.1 

SXDF850.3 

SXDF850.8 

SXDF850.il 

SXDF850.17 

SXDF850.119 



2.6 1 



3 +0.4 

, ,+0.3 
z,1 -0.1 
2.6 ' 



2.4 



+ 1.3 
0.1 
+ 0.4 
0.4 
> 2.0 
-1+0.1 
-0.1 



2 6 + 8 ° 
+10 

8 



27 



30 



+ 11 
-11 



46 
> 21 
2!) 



+ 29 



21 



+ 11 
-11 
+ 15 
- 14 



, 3+O.8 

2 3+ 1 ' 7 
2.4"*" ' 



„+3-9 
-2.0 
7+190.0 

> 0.9+ 23 „° 



7.7+ 



1.5 



-0.9 
+ 27.0 
-1.4 



2 2+ 4 ' 4 

4+1-5 
-0.6 

0. 



1.4"> 



+ 2.0 
-0.4 
0.9 



0.32 
> 2.2+f. 
0.6 4 



1-2.8 
-0.3 



400+ 130 
* uu -300 
400+300 
* uu -270 
410 + 680 
* lu -340 
, 3„„ + 32000 

+ 3900 
150 
+ 4700 
240 



> 160 
260 



X < 1 

x 2 « 1 
x 2 « 1 
x 2 « 1 
x 2 < 1 
x 2 « 1 



blackbody emission by a certain mass of dust) and ka can 
be extrapolated from an average Ki2s„„ = 2.64 m 2 kg -1 as- 
suming P = 1.5 ijDunne. Eales fc Edmundll2003h . 

Using Equation [21 the median dust mass implied by the 
850 ;Um observations is 9 x 10 8 Mq. Uncertainties in the dust 
masses are dominated by the uncertainty in the photometric 
redshifts (when these are used), Td and the flux density at 
850 /J,m. The uncertainty in k„ is a factor of a few and has 
not been included in the dust mass error bars, but the rel- 
ative dust masses in our sampl e will be correct i f the same 
value holds for all SMGs (see iBlain et all |2002| and refer- 
ences therein for a discussion) . Assuming that the maximum 
possible interstellar dust mass fo r a galaxy is about 1/50 
of its total baryonic mass (see lEdmunds fc Ealesl Il99ct ). 
the total baryonic mass of each SMG is estimated to be 
Mbary —5 x 10 11 Mq . For comparison, the baryonic mass of 
the Milky Way is about 10 11 Mq. 

The best-fitting SEDs are given in Fig.H while Table [5] 
gives the best-fitting value of Td for each SMG, as well as 
the derived Lfir, and SFR. 

Errors in the fit parameters (i.e. the position of the peak 



and the SED normalisation) are derived from the 68 per cent 
X 2 confidence interval around the best-fitting values. These 
errors have been carried through to all the other derived 
quantities in Table [5] Obviously, for galaxies with only two 
FIR photometric points, the SEDs are not well-constrained 
and the minimum reduced \ 2 values are significantly less 
than one. Photometric redshift errors have been folded into 
the uncertainties for derived quantities. We assume that our 
adopted model (a greybody with j3 = 1.5) is correct and so 
the quoted errors do not include a contribution from pos- 
sible systematic errors. Adopting a constant value of /3 is 
equivalent to assuming that the dust has similar properties 
in all galaxies. This is probably a reasonable assumption. As 
additional precise photometric data become available (from 
Herschel for example) /3 will be able to be constrained by 
direct fitting. Note, however, that changing the value of /3 by 
±0.5 (between physically plausible values) has the effect of 
changing the derived dust temperature by about ±5-10 K, 
which in turn affects the derived luminosities of the SMGs 
by ±20-40 per cent. 

The SMG population is therefore confirmed to be dom- 
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Figure 5. Apparent (observed) and restframe dust temperature 
distributions of an unbiased sample of SMGs. Each distribution 
has been calculated by adding a rectangle of area one centred 
on the source temperature (Ta or T d ), with a width given by 
the errors in the temperature (see Table [SJ to the histogram. 
The Ta distribution reveals that our survey selects SMGs with 
Ta ~ (8 ± 3) K. The T d distribution contains bumpy features 
around 12 K, 28 K, and 35 K and is found to be broader than we 
would expect by chance (see text). The underlying distribution 
is consistent with most SMGs having T d ~ 28 K but with a few 
being a bit hotter (~ 35 K). 

mated by massive (> 1O 12 M0, assuming Mtot /-Mbaryons — 
6), luminous (~ 2 x 10 12 L Q ) star-forming (SFR ~ 
400 Mq yr _1 ) galaxies with dust temperatures of ~ 35 K. In 
contrast, loca l starb urst galaxies with similar temperatures 
jDunne et alj l2000l : T d = 36 K) are not usually classified 
as ULIRGs, since they are about ten times less luminous 
(£fir ~ 10 11 L©). In addition, SMGs have more than ten 
times more dust than local starburst galaxies, which could 
be explained by any of the following reasons: (1) SMGs are 
hosted in galaxies ten times more massive than local star- 
bursts; (2) SMGs are more gas- and dust- rich for a given 
baryonic mass, being at an earlier stage in evolution with a 
higher gas fraction; or (3) the dust properties evolve with 
redshift and high-redshift dust has a higher emissivity and 
so produces more submm emission for a given mass of dust. 



5 DISCUSSION 

5.1 Trends in the data and consideration of 
selection effects 

We first consider selection effects in our data before looking 
for trends in the derived properties. It is not surprising that 
the S350 /Sgso colour appears to be independent of redshift 
(see Fig. [6}, owing to the small dynamic range in flux densi- 
ties of the 850 ^im-selected sources. The colour could evolve 



Figure 6. Observed S350 / $850 colour versus redshift using the 
values from Table [4] See inset for a description of the plot sym- 
bols. Sources with only redshift lower limits are omitted from 
this plot. The colour appears to be independent of redshift in 
these data (the median colour is 3.0), which is not surprising 
given the lack of dynamic range in the flux densities of the 
detected sources at 850 fim. Given a higher dynamic range in 
flux in future studies, we could attempt to trace SED evolu- 
tion with redshift by distinguishing between models. We have 
plotted some plausible models including (see inset fo r legend) an 
isoth ermal greybody with T d = 40 K, j3 = 1.5 (e.g. iBlain et al.l 
|2002|) ; an isothermal greybody T d = 34 K, (3 = 1.5 (this work); 
a 2-component typical IR-bright starburst galax y with T d = 20 
fc 40 K, = 2, and Af cold /A/ Warm = 30 (e.g. iPunne fc Ealesl 
1200 ll) ; a theoretical warm ULIRG T d = 24 & 50 K, /3 = 2, and 
M cold /M W arm = 5; a normal spiral (e.g. Milky Way, NGC891, or 
M51) with T d = 15 & 30 K, /3 = 2, and M cold /M warm = 100; 
a cold ULIRG (e.g. Arp220) with T d = 18 & 48 K, /3 = 2, and 
Afcold/Mwarm = 42 . 



with redshift, given a wider dynamic range. With this in 
mind, the relationship between dust temperature and red- 
shift is examined in Fig. [7] Notice that the spectroscopic 
and photometric redshift samples have about the same scat- 
ter and that the data deviate about the horizontal line in 
Fig. [7] (the sample median) by about ±5 K between redshifts 
of about 1 and 4 (except for a few outliers). It appears that 
the SMGs span a range of dust temperatures, with an ap- 
parent trend of T d increasing with z. A Spearman rank test 
reveals a strong correlation in the data with spectroscopic 
redshifts and a weaker correlation when the photometric red- 
shift sample is included. Is this due to selection effects or is 
this evidence for evolution of SMGs/ULIRGs with redshift? 

In order to test for selection effects in our data, we pop- 
ulate the Td-z parameter space of Fig. [7] by assuming rea- 
sonable distributions of SMG properties and apply 850 and 
350 /im cuts to the data to see if the trends seen in Fig. [7]can 
be reproduced by selection effects alone. The dust tempera- 
ture distribution of SMGs is assumed to be a simple Gaus- 
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Figure 7. Derived dust temperature versus redshift for the sam- 
ple of SMGs. See inset for a description of the plot symbols. The 
horizontal line represents the median dust temperature of the 
sample. The diagonal solid line indicates an observed dust tem- 
perature of 10 K, i.e. T d = 10 (1 + z)K, and appe ars to fit the 
data b etter; a similar conclusion was also noted bv lKovacs et al] 
(2006). This is not surprising, given that we find a constant 
•S l 35o/585o with redshift, as seen in Fig. [6] Combined with Monte 
Carlo simulations of the sample selection effects (see text), it ap- 
pears that the absence of 30 K sources seen at z > 3 is mainly 
due to a selection effect from SHARC-II, indicated by the shaded 
region. The lack of objects in the upper left is the result of evo- 
lution and shows that the 850 )im survey detects fewer warm 
SMGs in the low-redsh ift Universe than at high-rcdshift. See also 
IChapman et al.l ll2005h . 



sian with a me an of 35 K and a = 6K, consistent wit h ob- 
servations (e.g. iKovacs et al.|[2o5r3 ; iHuvnh et alj|2007t ). We 
create a /3 = 1.5 greybody SED with a temperature ran- 
domly selected from the temperature distribution, with a 
normalisation constrained so that Lfir is within the range 
of what we know about SMGs, i.e. a simple Gaussian distri- 
buti on with a mean of 6.7 x 10 12 L and a = 3.0 x 10 12 L Q 
(e.g. iPope et al.l [2005h . A redshift for each source is ran- 
domly selected from a simple uniform redshift distribution 
spanning < z < 5 to allow observed 850 and 350 fim 
fluxes to be calculated. We reject an SMG from our surveyed 
sample if its 850 and 350 /im flux densities are less than 
4 mjy and 25 mjy, respectively. This procedure is repeated 
100 times and the shaded region in Fig. [7] indicates the re- 
gion where selection effects preclude the detection of SMGs 
in our 350 /im survey of 850 ^m-selected SMGs. We miss 
sources in the lower right-hand region (high-redshift, low 
temperature) of Fig. [7] because they are not bright enough 
at 850 to make it into our SHADES sample. In addition 
some of the 850 /im-selected sources are too faint at 350 fim. 
On the other hand, there seems to be a dearth of sources in 
the data at low redshifts with high dust temperatures (and 
consequently, high luminosities), where we are not prone to 



Figure 8. Bolometric luminosity versus redshift for the sample 
of SMGs. The plot symbols are the same as in Fig. [6] except that 
arrows now indicate points with redshift lower limits. The shaded 
region indicates selection effects from SHARC-II. There appears 
to be a dearth of higher luminosity sources at low redshift: a hint 
of a strong luminosity evolution for SMGs with redshift. 



any selection effects (Fig. . The lack of sources in the up- 
per left of Fig. [7] is probably a consequence of the long- 
established strong lu minosity evolution of the su bmm pop- 
ulation (see fig. 2 in lSmail. Ivison fc Blainlll997l ). i.e. there 
are not as many sources at low redshift with temperatures 
> 35 K with sufficient luminosity to be detected by SHADES 
at 850 /im. Of course there could still be a systematic tem- 
perature dependence which might be seen by deep 450 and 
850 fim SCUBA-2 data through the detection of many more 
450 fim sources at fainter flux densities. 

This idea can be further investigated when the relation- 
ship between Lfir and redshift is examined (Fig. [HJ, since 
there is an intrinsic correlation between temperature and 
Lfir. As before, SHARC-II is insensitive to the shaded re- 
gion in Fig.|8](i.e. misses SMGs cooler than T d /(l+z) ~7K). 
Due to the increased volume at low redshifts, we would ex- 
pect to see about two high-luminosity sources between red- 
shifts of and 1. Thus there appears to be a hint of a dearth 
of higher luminosity SMGs at the lowest redshifts compared 
to higher redshifts that cannot be explained by selection or 
volume effects. 

These results are difficult to make sense of unless a 
strong SMG evolution with redshift is invoked. Millimetre- 
wave cameras such as MAMBO, Bolocam and AzTEC might 
be sensitive to cooler SEDs and could perhaps select SMGs 
to fill in this region of the (Lfir, z) plane. There also seems 
to be a lack of SMGs at low redshifts and high luminos- 
ity, indicating that SCUBA-selected SMGs are intrinsically 
more lumino us at high red s hifts. This is consistent with what 
was found bv llvison etafl l|2002l ) and lPope etafl (|2006l ) and 
suggests a strong evolution of SMGs/ULIRGs with redshift. 
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Thus, Fig. [7] is just a reflection of the known effect seen in 
Fig. [8] i.e. that the most luminous SMGs are at higher red- 
shifts. This means that they must be evolving very steeply 
with redshif t (i.e. 'down-sizing'). This i dea has already been 
explored in IWall. Pope fc Scottl (|2007h . Larger samples of 
SMGs with robust redshifts and a wider dynamic range in 
350 and 850 /im flux densities and luminosity are needed to 
test this idea further and to fully describe the evolution of 
the luminosity function of SMGs. 
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6 CONCLUSIONS 

We have performed follow-up mapping at 350 /im us- 
ing SHARC-II of several 850 /mi-selected sources from 
SHADES. In total there are flux density measurements for 
about 25 per cent of the SHADES SMGs, of which this work 
has provided 21 new flux density constraints. 

The combination of 350 and 850 /im, 1.1, and/or 1.2 mm 
photometry and spectroscopic or photometric redshifts have 
provided estimates of Td and Lfir for each 350 /im-observed 
SHADES SMG by fitting a modified blackbody function to 
the data. The SMGs population is confirmed to be dom- 
inated by very dusty (Md us t - 9 x 10 s M0), luminous 
(L Fm ~ 2 x 1O 12 L ) star-forming (SFR ~ 400M o yr~ 1 ) 
galaxies with dust temperatures of ~ 35 K. We have mea- 
sured the temperature distribution of SMGs and find that 
the underlying distribution is slightly broader than implied 
by the error bars, and that most SMGs are at 28 K with a 
few hotter. 

With the data in hand, it is very difficult to understand 
whether we are just seeing a combination of selection effects 
and the known luminosity evolution of SMGs, or whether 
the SEDs might also systematically change with redshift. 
Larger samples of SMGs with a broader selection function 
will enable this to be fully investigated, e.g. with a deep 
450 /mi survey, a primary goal of the upcoming SCUBA-2 
Cosmological Legacy SurvejQ. 

Using these data, we have also been able to estimate 
the 350 /im source counts at 200 deg -2 < N(> 25mJy) < 
500 deg -2 . Large blank-field surveys at 350 and/or 450 /mi 
will allow us to study a statistically significant population 
of lower redshift and/or hotter dust temperature galaxies, 
responsible for emission nearer to the peak of the extra- 
galactic background than those selected at 850 /im. These 
studies should provide a link between galaxies detected at 
~ 1 mm and those selected in the mid-IR with Spitzer. Such 
comparisons should allow us to investigate the evolution of 
FIR emitting galaxies and how they contribute to the back- 
ground at different wavelengths. 
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APPENDIX A: EFFECTIVE EXPOSURE TIMES 

When describing astronomical data, it is common and 
straightforward to give total observation times, or perhaps 
in chopped systems the total time on source. These times 
give an indication of the amount of effort expended to col- 
lect the data. However, because the atmosphere is far from 
transparent at submm wavelengths (even in the atmospheric 
windows) and the transmission is so weather dependent, the 
total exposure time does not give a good indication of data 
quality or expected noise level in the final output. It is more 
useful to refer to the effective total exposure which is the 
noise-squared weighted sum of exposure times referred to a 
completely transparent atmosphere. 

The optical depth of the atmosphere in a given direction 
is determined from the wavelength dependent optical depth 
at zenith, t\, and the cosecant of the zenith angle, 6, of 
observation, called the airmass, A. Thus, the strength of an 
observed astronomical flux density, So is: 

So = S e~ rMAW . (Al) 

At 850 (im, the zenith optical depth, rsso, ranges from, 
say, 0.16 in very good weather to 0.7 in bad weather, and the 
optical depth at 350 or 450 /im is typically 7 times higher. 
The first step in combining data collected at different times 



of night with different values of t\ and A is to correct 
So for atmospheric ab sorption by multiplying by e TA A (see 
I Archibald et al.ll2002l ). 

After correcting for atmospheric transients, etc., but not 
opacity, the output noise at the detector is fairly independent 
of atmospheric opacity, at least for SCUBA operating on the 
JCMT on Mauna Kea at 450 and 850 /Ltm, and for SHARC- 
II operating at 350 ^m at the CSO also on Mauna Kea. 
Therefore, the effective noise of an observation is inflated 
exponentially in t\ and A. The sum 

t cff = j2 At ° e2Tx(t)A ( A2 ) 

is the integration time that would be required with a com- 
pletely transparent atmosphere to get the same noise level 
as a given noise-squared weighted sum of observations, each 
of size A t - We call this the effective exposure time and find 
that tracking this quantity during observations is an effec- 
tive way to monitor experimental progress. 

In practice, the optical depth as a function of time is 
obtained from automated measurements made at the CSO 
at 225 GHz every 10 minutes. The optical depths at submm 
wavelengths are very well correlated with these values and 
one multiplies T225 GHz by about 4 and by 24 to get optical 
depths for SCUBA at 850 and 450 fj, m, respectively, and b y 
26 to get r for SHARC-II at 350 /im (| Archibald et alj|2002h . 

As a numerical example, observation of a source which 
rises through transit with an average airmass of 1.3 and 
t 225 GHz = 0.05 for 4 hours results in effective exposure 
times of 2.4 hours at 850 /jm, 20.7 minutes at 450 /xm with 
SCUBA, and 9.5 minutes with SHARC-II. The shorter effec- 
tive times at shorter wavelengths are a consequence of the 
fact that the atmosphere is not very transparent at those 
wavelengths. This is why space experiments like Herschel- 
SPIRE and BLAST are useful even with comparatively 
small apertures, and why neither of those experiments car- 
ries an 850 /im channel. Experience has shown that effective 
exposures of between 500 and 600 seconds with SHARC-II 
are often sufficient to detect faint extragalactic high-redshift 
sources. 
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